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Abstract 

We consider elastic quark-quark scattering at high energy and fixed transferred mo- 
mentum. Performing factorization of soft gluon exchanges into Wilson line expectation 
value we find that there is one-to-one correspondence between high energy asymptotics 
in QCD and renormalization properties of the so called cross singularities of Wilson lines. 
Using this relation we show that the asymptotic behavior of the quark-quark scattering 
amplitude is controlled by a 2 x 2 matrix of the cross anomalous dimensions. We evaluate 
the matrix of cross anomalous dimension to two-loop order and study the properties of 
the obtained expressions to higher loop order. 
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1. Introduction 



Recent experimental data from Tevatron and HERA renewed our interest to theoretical in- 
terpretation of the well known phenomena in the physics of strong interactions like growth 
of the total hadronic cross section at high energy and increasing of the structure function of 
deep inelastic scattering at small values of the Bjorken variable x. In both cases we deal with 
the asymptotic behavior of hadronic scattering amplitudes at high energy s and fixed trans- 
ferred momenta t the so called Regge behavior in high-energy QCD This problem is one 
of the longstanding in QCD and previous attempts to solve it have led to the discovery of the 
BFKL pomeron ||. However, found in the leading logarithmic approximation, cr s C 1 and 
a s logs ~ 1, the BFKL pomeron leads to the expressions for the scattering amplitudes which 
violate the unitarity Froissart bound. To preserve the unitarity of the S— matrix, one has to 
go beyond the leading logarithmic approximation calculating the scattering amplitudes. The 
problem turns out to be very complicated and the following approaches have been proposed: 

• Generalized leading logarithmic approximation |3|, |J], in which we start with the lead- 
ing logarithmic result and add "minimal" number of nonleading corrections in order to 
unitarize the result; 

• Two-dimensional reggeon field theory |J, in which the high-energy scattering of gluons 
is described by the S— matrix of the effective two-dimensional field theory corresponding 
to dynamics of transverse gluonic degrees of freedom; 

• Quasielastic unitarity approximation ||, in which one takes into account only contribu- 
tions of soft gluons with transverse momenta much smaller hadronic scale. 

These schemes are approximate and they were developed to resum special class of corrections 
which are expected to determine the high-energy asymptotics of the scattering amplitudes. In 
each case, QCD is replaced by simplified effective theory which inherits a rich structure of the 
original theory. 

Recently 17], 0, a progress has been made using the generalized leading logarithmic approxi- 
mation. It was found that in this approximation high energy QCD turns out to be a completely 
integrable model equivalent to the one-dimensional Heisenberg magnet of spin s = 0. This opens 
the possibility to apply powefull quantum inverse scattering method and calculate the hadronic 
scattering amplitudes by means of the Generalized Bethe Ansatz ||. 

There were different proposals for the effective two-dimensional reggeon field theory |5|, [| 
and the hope is that the model can be solved exactly. Although to the lowest order of PT the 
effective reggeon theory reproduces the QCD expressions for the scattering amplitudes, the all 
order perturbative solution is not available yet. 

The present paper is devoted to the calculation of the hadronic scattering amplitudes in the 
quasielastic unitarity approximation 0. In this approximation, we use the parton model and 
approximate the hadron-hadron interaction as the scattering of partons. Partons interact with 
each other by exchanging soft gluons with energy which is much smaller than parton energy, s, 
but much bigger than the QCD scale Aq C d in order for the perturbative QCD to be applicable. 
Being calculated to the lowest orders of PT, the parton-parton scattering amplitudes get large 
perturbative corrections (a s logs log t) n due to soft gluon exchanges || and these Sudakov 
logarithms need to be resummed to all orders of PT. The resummation of leading logarithmic 



corrections can be performed using the "evolution equation" technique p, |10[ but we don't 
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have a regular way to resum and control nonleading logarithmic corrections to the scattering 
amplitude. A new approach which enables us to take into account systematically both leading 
and nonleading logarithmic corrections to the parton-parton scattering amplitudes has been 



proposed in 11 1. In this approach, the partonic scattering in QCD is described by the dynamics 



of the gauge invariant collective variables defined as 

W = W{d,C 2 , ...,C n ) = (0|T tr Rl Pe ig ^ dx ' A(x) ■ ■ ■ tr^Pe^^^O) , (1.1) 

the so called Wilson loop functions. Here, the path-ordered exponentials are evaluated along 
the closed paths C\ . . . C n in Minkowski space-time and the gauge fields are defined in the 
irreducible representations R\ . . . R n of the SU (N) gauge group. One should stress that it is not 
for the first time when one dealt with this object. There were attempts [12 in 80's to study the 



nonperturbative dynamics of QCD in terms of Wilson loops. The interest to the Wilson loops 
was renewed after it was recognized fl3| that the same loop functions enter into consideration 



in perturbative QCD when we study the effects of soft gluons in hardonic processes. Interacting 
with soft gluons described by gauge field A^x), partons behave as relativistic charged classical 
particles and the path-ordered exponentials appear as eikonal phases of partons. The integration 
paths C\ . . . C n have a meaning of semi-classical trajectories of partons and the choice of the 
representations R\ . . . R n is fixed by their charges. In particular, R is fundamental and adjoint 
representation of the SU(N) group for quarks and gluons, respectively. 

The paper is organized as follows. In sect. 2 we consider the amplitude of near forward 
elastic quark-quark scattering and show its relation with the renormalization properties of the 
so called cross singularities of Wilson loops. In sects. 3 and 4 we analyze the cross singularities 
of Wilson loops to the lowest orders of PT and find two-loop expression for the matrix of cross 
anomalous dimensions. The properties of the obtained expressions which can be generalized to 
all orders of PT are studied in sect. 5. The details of the calculations of the two-loop Feynman 
diagrams contributing to the Wilson lines are given in Appendices A and B. Large N limit of 
the Wilson line is analysed in Appendix C. 



2. High energy scattering in QCD 

Let us briefly summarize the results of the paper |TI| concerning the relation between high 
energy behavior of the scattering amplitudes in QCD and the renormalization properties of 
the Wilson loops. We consider the near forward elastic quark-quark scattering in the following 
kinematics: 

s, m 2 > -t > A 2 > A 2 CD . 

Here, s = {pi + P2) 2 is the invariant energy of quarks with mass m, t = (pi — p[) 2 is the 
transferred momentum and A 2 is the IR cutoff (characteristic hadronic scale) which is introduced 
to regularize infrared divergences. Notice that we don't need to fix the relation between s and 
m and quarks may be heavy. The scattering amplitude T-j J contains the color indices of 
the incoming and outgoing quarks. The quarks interact with each other by exchanging soft 
gluons with the total momentum q — pi — p' x . In the center of mass frame of incoming quarks, 
for s ^> —t, the transferred momentum has the light-cone components q = {q + ,q~,q) with 
q ± = ^(g° ± q 3 ) = and q 2 = —q 2 = t. Since the transferred momentum is much smaller 
than the energies of incoming quarks, the only effect of their interaction with soft gluons is 
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the appearance of additional phase in the wave functions of the incoming quarks. This phase, 
the so called eikonal phase, is equal to a Wilson line Pe 1 Ic dx ' A ^ defined in the fundamental 
representation of the SU (N) group and evaluated along the classical trajectory C of the quark. 
We combine the eikonal phases of both quarks and obtain the representation for the quark-quark 
scattering amplitude as [ITT] 




sinh 7 / d l z e- lz '- q ' W?f (7, z 2 X 2 ) , t = -q 2 , (2.2 



where the line function W^ 3 is given by 



( 7 ,^A 2 ) = (0|T 



Pe 



in 



J°° davi-A(via) 



Pe 



ig df3v2-A(v2p+z) 



|0>. (2.3) 
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Here, two Wilson lines are evaluated along infinite lines in the direction of the quark velocities 
v i = Vi/ m an d i>2 — Vil m and the integration paths are separated by the impact vector 
z = (0 + , 0~, z) in the transverse direction, v\ ■ z — v 2 • z — 0. 

Being expanded in powers of the gauge field, expression ( |2.2j ) reproduces the quark-quark 
scattering amplitude in the eikonal approximation. Integration over two-dimensional impact 
vector z was introduced in ( |2.2| ) to ensure the total transferred momentum in the t— channel 
to be q = (0 + ,0~,q). Indeed, it follows from (|2.2|) that the total momentum of gluons k to t. 
propagating in the t— channel is restricted by the following 5— functions: 

S(k to t. ■ v^Siktot. ■ v 2 )S(k to t. - q) = . 7 6 (k^ )S(ktot )S(k to t. ~ q) = . ! ^ 4 \k tot , - q) , 

sinh 7 ' ' sinh 7 



where the first 6— function comes from integration over a in ( |2.3|) , the second one from inte- 
gration over j3 and the last one from integration over z in ( |2.2|J . To compensate the additional 
factor, sinh 7, in the r.h.s. of this relation, the same factor was introduced in the expression 

The scattering amplitude (|2.2|) depends on the quark velocities Vi and v 2 , the transferred 
momentum q and IR cutoff A. These variables give rise to only two scalar dimensionless 
invariants: (17^2) = s/m 2 and t/X 2 as explicitly indicated in ( |2.2j ). The s— dependence of the 
amplitude comes through the dependence on the angle 7 between quark 4— velocities V\ and v 2 
defined in Minkowski space-time as 

(v\ • v 2 ) — cosh 7 . (2.4) 

For arbitrary velocities V\ and v 2 , the angle 7 may take real and complex values and in the 
limit of high energy quark-quark scattering (s ^> m 2 ) we have 

7 = log^»l. (2.5) 

To find the asymptotic behavior of the scattering amplitude ( |2.2|) one needs to know the vacuum 
expectation value of the Wilson lines ( f2.3| ). Let us consider first a special case of the elastic 
electron-electron scattering in QED. 
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2.1. High-energy scattering in QED 



One may use the expression ( |2.2|) for the scattering amplitude in this case, but the calculation 
of the Wilson line ( j2.3[ ) is simplified because for abelian gauge group a path-ordered exponent 
coincides with an ordinary exponent. Moreover, since photons don't interact with each other, 
the calculation of Wilson lines in QED is reduced to the integration of a free photon propagator 
D^uix) along the path and the result of calculation is 



Wi 



QED l 



(7,^ 2 A 2 )=exp [^Ljdx* Jdy v D, v {x-y) ) = (\ 2 z 2 /A)^ coth \ (2.6) 



where A is a fictitious photon mass. After substitution of this relation into ( |2.2| ) and integration 
over impact vector we get the expression for the scattering amplitude in QED H 

s t\ . 2 cosh7 /-A-^ coth 7 r(l + igcoth7) 



M^J^^UJ ra-igcth,; - (2 ' 7) 

In the limit s/m 2 — > oo or, equivalently 7 — > 00, the expectation value ( j2.6|) does not depend 
on s and, as a consequence, the scattering amplitude has the asymptotics T/T Born ~ s° with 
^Born = the amplitude in the Born approximation. This implies that photon is not 

reggeized in QED Hl4fl . 

2.2. Wilson lines in perturbative QCD 

The evaluation of the Wilson line in perturbative QCD is much more complicated than in QED. 



To explain the approach |TT| for the evaluation of the Wilson line ( |2.3|) , let us consider as an 
example the one-loop calculation of W in the Feynman gauge. Using the definition ( |2.3| ), we 

get 



W^ loop = /®/+(r®r)-^-r( J D/2-l)A 4 - D da dp 

47T ' J— 00 J — 00 



{v 1 a-v 2 P) 2 + z 2 + iO} D / 2 - 1 ' 

(2.8) 

where a direct product of the gauge generators defined in the fundamental representation takes 
care of the color indices of the quarks. In this expression, gluon is attached to both Wilson 
lines at points V\a and v 2 P + z and we integrate the gluon propagator ViV^D^iy x a — v%ft — z) 
over positions of these points. To regularize IR divergences we introduced the dimensional 
regularization with D = 4 — 2e, (e < 0) and A being the IR renormalization parameter. There 
is another contribution to W^i-ioop corresponding to the case when gluon is attached by both 
ends to one of the Wilson lines. A careful treatment shows that this contribution vanishes. 
Performing integration over parameters a and (3 in (|2.8|) we get 



Wi-ioo P (7,A 2 ? 2 ) = J®/ + {t a ®t a ) — (-mcoth7)r(-£)(7rA 2 f 2 ) £ . (2.9) 

71 

The integral over a and f3 in ( |2.8|) has an infrared divergence coming from large a and (5. In 
the dimensional regularization, this divergence appears in W^i-ioop as a pole in (D — 4) with 
the renormalization parameter A having a sense of an IR cutoff. After renormalizaion in the 
MS— scheme, this result coincides up to color factor with the lowest term in the expansion of 
the analogous expression (|2.6|) in QED in powers of e 2 . 
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The evaluation of the Wilson line (|2.3| ) in QCD is based on the following observation. The 
one-loop expression (|2.9|) for the line function is divergent for z = and e < 0. This divergence 
has an ultraviolet origin because it comes from integration over small a and (3 in ( |2.8| ), that is 
from gluons propagating at short distances (via — v 2 /3) between quark trajectories as a, (3 — > 0. 
One should notice that the ultraviolet (UV) divergences of Wilson lines at z = have nothing 
to do with the "conventional" ultraviolet singularities in QCD. For z = 0, the integration 
paths of Wilson lines in (|2.3|) cross each other at point and UV divergences, the so called 



cross divergences, appear when one integrates gluon propagators along integration path at the 



vicinity of the cross point |TJJ . We note that the cross divergences do not appear in Wi_i oop for 
nonzero z because, as it enters into the integral (|2.8|) , nonzero z 2 regularizes gluon propagator 
at short distances as a, (3 — > 0. 

Thus, z 2 has a meaning of an UV cutoff for the one-loop line function W and the depen- 
dence of the original line function W(^, \ 2 z 2 ) on z 2 is in one-to-one correspondence with the 
dependence of the same Wilson line but with z = on the UV cutoff which we have to introduce 
to regularize the cross singularities. This property suggests the following way for calculation of 
the z— dependence of the Wilson lines (JO|). First, we put z = in the definition ( |2.3| ) of the 
line function and introduce the regularization of cross singularities of W. Then, we renormalize 
cross singularities and identify UV cutoff with impact vector as /i 2 = 1/z 2 . We conclude that 
the z— dependence of the line function W and, as consequence, the t— asymptotics of the scat- 
tering amplitude (|2.2|) are controlled by the renormalization properties of the cross singularities 
of Wilson loops. 

2.3. Renormalization properties of Wilson lines 

The renormalization properties of the loop functions in QCD have been studied in detail |fT5| 
and can be summarized as follows. Being expanded in powers of gauge potential, W{C\, . . . , C n ) 
can be expressed as multiple integrals of the Green functions along the paths C\ ... C n . The 
loop function is not a well defined object because it contains different kinds of divergences. 
First of all, it has the conventional ultraviolet (UV) divergences of the Green functions in 
perturbative QCD. Moreover, as it was first found in |Tj|, the Wilson loops have very specific 
UV divergences which depend on the properties of the integration paths C in the definition 
( |2.3| ). For instance, the expansion of the loop function to the lowest order of PT involves 
the contour integral § c dx^ § c dy u D^ v {x — y) which is potentially divergent at x = y due to 
singular behavior of the gluon propagator at short distances. The divergences come either 
when y approaches x along the path C or when the path crosses itself and both x and y 
approach the cross point. In both cases divergences have an UV origin because they come 
from integration at short distances. In the first case, UV divergences can be absorbed into 
renormalization of the coupling constant and gauge fields provided that the integration path is 
smooth everywhere. Otherwise, if the path is smooth everywhere except of a one point where 
it has a cusp, additional UV divergences, the so-called "cusp" singularities, appear which can 



be renormalized multiplicatively [17]. The corresponding anomalous dimension, r cusp (7,g), the 
cusp anomalous dimension is a gauge invariant function of the coupling constant and the cusp 
angle. If the integration path crosses itself (or the paths C\ . . . C n cross each other), additional 



"crossed" singularities appear in the loop function [15[ . 



It was found JTB) that these renormalization properties of the Wilson loops which look like 



their undesired features are of the most importance in perturbative QCD. The same function, 
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the cusp anomalous dimension of a Wilson line, appears when one studies the infrared (IR) 
singularities of on-shell form factors or the velocity-dependent anomalous dimension in the 
effective heavy quark field theory. The anomalous dimensions of the Wilson lines govern the 
IR logarithms in the same manner as the anomalous dimensions of the local composite twist- 2 
operators control the collinear singularities one encounters studying the structure functions 



of deep inelastic scattering. Moreover, there is a remarkable relation JL1]] between the high- 
energy behavior of the scattering amplitudes in the quasielastic unitarity approximation and 
the renormalization properties of the "cross" divergences of the loop function W. 

2.4. Cross anomalous dimension 

The renormalization properties of the cross divergences were studied in |15|. It was shown 



that the loop function (|1.1| ) with all Wilson loops defined in the fundamental representation 
of the SU(N) group (Ri — ■ ■ ■ — R n — F) and paths having a cross point, is mixed under 
renormalization with loop functions evaluated along the integration paths which are different 
from original paths at the vicinity of the cross point. As an important example we consider 
the Wilson lines (Wx) l J and (W^/ defined in fig.l. Although the integration paths in fig.l 
have an infinite length and they are open we could get two Wilson loops by projecting the line 
functions (Wi) % J and {W^fJ onto the color structure S^jSju. 

We choose the Wilson lines in fig. 1(a) and 1(b) to be defined in the fundamental represen- 
tation of the SU(N) gauge group. According to the general analysis []15| , the line function of 



fig. 1(a) has cross singularities and under their renormalization it is mixed with the Wilson line 
of fig. 1(b). As a consequence, the renormalized line functions W\ and W2 satisfy the following 



renormalization group (RG) equation 



{^ + (3{g) ^) Wa = - r "°»( 7 ' 9)Wb ' a ' 6=1 ' 2 ' (2 ' 10) 

where W a = (W a )^ and jj, is renormalization point. Here, r cross (7,g) is the cross anomalous 
dimension which is a gauge invariant 2x2 matrix depending only on the coupling constant and 
the angle 7 between lines at the cross point. It is important to notice that the RG equation 
( 2.10| ) is valid only in the fundamental representation of the gauge group and its generalization 



to another representations is uknown. 



2.5. Evaluation of the scattering amplitude 



We recall that the evaluation of the scattering amplitude ( |2.2[ ) is based on the following property 
of the line function W. Using the definition (|2.3| ) of W we introduce into consideration the 
Wilson line W± in fig. 1(a). The only difference between them is that the integration paths in 
( |2.3| ) do not cross each other for nonzero z. For z 7^ the Wilson line W(7, i^A 2 ) is UV finite 
but for z = it coincides with W\ and has cross divergences. The distance between quark 
trajectories, z 2 , in the definition (|2.3|) of W plays a role of UV cutoff for cross singularities 



and the z— dependence of W is the same as the dependence of the Wilson line W\(j, /x 2 /A 2 ) on 
the renormalization parameter fi. Solving the RG equation ( |2.10| ) for ^1(7, fi 2 /X 2 ) with the 
boundary conditions, Wi(j, 1) = SinSfj and W^(7, 1) = we substitute fi 2 = 1/z 2 and 
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find the following expression for the scattering amplitude [|TTJ 



r if (^2> % J = sinh 7 / (A n (7, ^A 2 ) <W<^' + A 12 ( 1} ?\ 2 ) <^y) , (2.11) 

where An and v4 12 are elements of the 2x2 matrix 

A qq ( 7 ,z 2 A 2 ) = T exp ( - ^ — r cross ( 7 ^(r))j. (2.12) 

Notice, that the matrices r cross (7, g(r)) do not commute with each other for different values 
of r and this makes difficult to evaluate the time-ordered exponent in fl2.12| ). The expression 



(|2.11| ) for the quark-quark scattering amplitude can be easily generalized to describe the quark- 
antiquark scattering. We notice that antiquark with a 4-velocity v 2 can be treated as a quark 
moving backward in time with velocity — v%. As a result, to get the quark- antiquark scattering 
amplitude one has to replace one of the quark velocities v 2 — > —v 2 in expressions (j2.11|) and 
( |2.12| ). Since the scattering amplitude (|2.11|) depends on the angle 7 between velocities, this 
transformation corresponds to the replacement 7 — > vk — 7 in the expression ( |2.12| ) for the 
matrix A, or equivalently in cross anomalous dimension: 

Aqq(7, z 2 ^ 2 ) = T exp ^- jf ^r cross (MT - 7, g{f))\ , 

where 7 is the angle between velocities of incoming quark and antiquark. 



The expression (2.11) takes into account all logt and logs corrections to the scattering 



amplitude. We conclude from ( |2.12j ) that it is the matrix of the cross anomalous dimensions 



which governs the high energy behavior of the scattering amplitude. As it follows from (|2.11|) 



the s— dependence of the amplitude comes from the 7— dependence of r cross (7, g) while the 
t— dependence originates from the evolution of the coupling constant. To find the asymptotic 
behavior of the scattering amplitude (|2.11|) we need to know the large 7— behavior of r cross (7, g) 
to all orders of perturbation theory. To this end, we perform the one-loop calculation of 
r c ross(7)fl0 i n t ne next section and try to find the properties of this matrix which can be 
generalized to higher loop orders. 



3. One-loop calculation of the cross anomalous dimen- 
sion 

To find the cross anomalous dimension we evaluate the line functions of fig. 1(a) and 1(b) to 
the lowest order of PT and substitute them into ( |2.10|) . In the Born approximation the line 
functions have trivial expressions 

W[ 0) = 6 n 5 fj = \1), W 2 (0) = 5 ifj 6 fi = 1 2) , (3.13) 

where |1) and |2) denote two states in color space. To one-loop order the line functions are 
given by path integrals similar to that in (|2.8|) . To evaluate the integrals we introduce the 
following parameterization of the integration paths of fig.l, 

C\ = via , Cn = v 2 /3 , Cm = vi/3 , Ci V = v 2 a , (3.14) 
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where < a < oo and — oo < (3 < and where v\ and v 2 are quark velocities. 

The integration paths in fig. 1(a) and 1(b) cross each other at point and go to infinity along 
the vectors v% and t> 2 . As a consequence, the line functions W\ and W 2 have UV divergences 
due to the cross point and IR divergences due to the infinite length of the paths. To regularize 
the cross singularities we use the dimensional regularization with D = 4 — 2s and in order to 
regularize IR divergences we introduce a fictitious gluon mass A. We fix the Feynman gauge 
and define the gluon propagator in the coordinate representation as D fiu (x) = g lxv D{x) with 

\ ■ f d D k e~ ikx 1 /•-*» D/2 _ 2 f 2 X 2 \ 



(27r) D A; 2 -A 2 + iO An D / 2 

After the decomposition of the integration path of fig.l( a) as C — C\ -\- Cxi ~l~ Cm Civ we 
get the set of Feynman diagrams of fig. 2 which contribute to one- loop corrections to the line 
function W\. Then, the contribution of the diagram of fig. 2(1) to W\ is given by 

(ig) 2 t"it a fj / (hp / dy„ D^(x - y) = -g 2 ^^)^ ®t a da df3 D( Vl a - v 2 @) , 

where integration is performed over positions of gluons on the paths C\ and C\\ and the direct 
product of the gauge group generators takes care of color indices of the line function. Substi- 
tuting the gluon propagator (ft.l5Q and integrating over a and (3 we get the expression for the 



"Feynman integral" 

F 1 = {igf f^dx, f c Jy v D^(x-y) = ^^ 7 coth 7 . (3.16) 

Here, the cross singularity comes from the integration over small x — y and it appears in the 
final expression as a pole in e. Calculation of the remaining diagrams of fig. 2 is analogous and 
gives the following results 

a s /47r/i 2 yr(l+£) . _ a s (An^\ £ T(l + e) 

V = --[-#-) ^^(--7)coth 7 , F S = -2F 4 = --^J (3.17) 

We notice that there are simple relations between these expressions: -^2(7) = —Fiiin — 7), 
F3 = -Fi(7 = 0). In Appendix A we explain their origin and show that analogous relations exist 
between two-loop Feynman integrals. 

To evaluate the one-loop correction to the line function W\ we take the expressions ( |3.16| ) 



and ( 3.17 ) and sum them with appropriate color and combinatoric factors: 



M l) = tUh( 2F ± + 2F ?) + Cf <\;<\r,C2F, + 4F 4 ) , 

where Cp = t a t a is the quadratic Casimir operator in the fundamental representation of 
the SU(N) gauge group. This expression contains a pole in e which is subtracted in the 
MS— scheme. Finally, the renormalized one-loop expression for the line function W\ is given by 

Wf> = fog z7rcoth T . (3.18) 

The calculation of the line function W2 involves the integration over the path of fig. 1(b) which 
is different from the path of fig. 1(a) only at the cross point. The contribution of the individual 
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Feynman diagram of fig. 2 to the line function W\ has the form of the product, W — F ■ C • w, 
of the corresponding Feynman integral (F), the color factor (C) and the combinatorial weight 
of the diagram (w). The difference between the paths of fig. 1(a) and 1(b) at the cross point 
is not important when one calculates the Feynman integral F and the combinatorial weight 
w but it does affect the color factors of the corresponding diagrams. This means, that to any 
order of perturbation theory the same Feynman integrals contribute to both line functions and 
the only difference between W\ and W 2 is that integrals enter with different color factors. In 
particular, using the expressions (|3.16|) and ( p. 17] ) we get the one- loop expression for the line 
function W 2 as 



W 2 {1) = C F Svfo&Fi + 4F 4 ) + ty%(2F 2 + 2F- 



3 j 



After substitution of the Feynman integrals ( |3.16|) and (|3.17|) into this expression we subtract 
a pole in the MS scheme and obtain the renormalized line function 



W 2 (1, = —\og^ -Cf (5i/, -^( 7 coth 7- 1) +4,^(7 coth 7- m coth 7- 1) . (3.19) 



The obtained one-loop expressions (|3.18| ) and ( |3.19|) for the line functions are valid in an 



arbitrary representation of the SU(N) group after proper redefinition of the Casimir operator 
Cp. The reason why the fundamental representation plays a special role comes from the fact 
that the direct product of the gauge group generators can be decomposed into the sum of 
invariant tensors 

Vki = -^M + \^jk- (3.20) 
After substitution of this decomposition into ( pj.18 ) and ((57E3) we find that one-loop correction 



to W\ contains a color structure of W 2 in the Born approximation, eq.( |3.13j ), and vice versa. 
As a consequence, both line functions are mixed with each other under renormalization. In an 
arbitrary representation of the gauge group we do not have a relation similar to ( |3.20| ) and the 
renormalization properties of the line functions become complicated. 

Substituting the one-loop results, ( |3.18| ) and ( |3.19| ), for the line functions into the RG 
equation (|2.10|) we find the one-loop expression for the matrix of cross anomalous dimension: 



r ^(^g) = — r( 7 ) , (3.21) 

p/ \ _ I — ]f coth 7 Z7rcoth7 \ 

~ ^ -7coth7 + 1 + nrcoth7 A^(7coth7 - 1) - f coth7 )' 

Here, 7 is the angle between quark velocities v± and v 2 in Minkowski space-time defined in 
(12! 



3.1. Properties of the one-loop r cross 

The obtained expression (3.21) for one-loop T CTOSS obeys the following interesting properties. 
All elements of the matrix T CTOSS have nontrivial imaginary parts which change their sign under 
the replacement 7 — > iir — 7: 

Im r cross (z7r - 7, g) = -Im r cross (7, g) . 
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As was shown in sect. 2. 5, this transformation of the angle allows us to relate quark-quark and 
quark- ant iquark scattering amplitudes. Calculating the determinant of the matrix (3.21) we 
obtain the expression 

det r cross ( 7 , g) = (^) V coth 2 7 ( 1 - jL) , (3.22) 



71 



which, first, has zero imaginary part and, second, is invariant under 7 — > iix — 7. This property 
implies that the elements of the matrix (3.21) are "fine tunned". An additional argument in 
favour of this observation comes from the consideration of the large 7 behavior of r cross . In the 
limit of large 7 different elements of the matrix have the following behavior 

r 11 ~r 12 ~oh°) r 12 ~r 22 ~oh) 

cross cross ^ \ I 1 1 cross cross ^ \ 1 1 ) 

which leads to the asymptotics det T 

cross ~ ^(t) ^ general. However, the one-loop result 
( |3.22| ) implies that the leading term 0(7) vanishes and the determinant has the asymptotics 



det r cross ~ O(j ) while trr cross ~ The eigenvalues of the matrix T(7) satisfy the 

characteristic equation 

T| - r± f A^7 coth 7 - iV - — coth 7 J + (tt coth 7 ) 2 = . (3.23) 

Solving this equation we find that in the large 7 limit both eigenvalues have positive real part 
(Re r± > 0). Moreover, one of the eigenvalues of matrix T is much larger than the second one: 

r + = W7 _( w + £)_^^-. + ( 7 -), 

AT 2 — 1 

r- = vr 2 ^ 7 - 1 + 0(7- 2 ), ( 3 - 24 ) 
with 7 = log(s/m 2 ) in the limit s ^> m 2 . 



3.2. Asymptotic behavior of the scattering amplitudes 



Let us substitute the one-loop expression (3.21) for the matrix r cross into ( |2.1l ) and (2.12). 
One-loop matrices r cross (7, g{r)) commute with each other for different r and we can omit 
T— ordering in the definition ( |2.12|) of the matrix A. After diagonalization of the one-loop 
matrix T(7) by a proper unitary transformation we get 



2 2 

z 



r+-r\ 



cxp 



-r_ 



^dTa s {T)\ r_-r( 7 ) 
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r_-r. 



cxp 



1 / z dr otAr) 
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where the eigenvalues F± were defined in ( |3.23|) . Finally, we find the following expression for 
the scattering amplitude fl2.11|) 



5 .,. 5 .,. T (0) +t a it a j T (8) 



(3.25) 



where two invariant amplitudes correspond to the exchange in the t— channel by states with 
the quantum numbers of vacuum and gluon and are given by 



T ( ) 



smii7 r + r_ 
t - rv -r_ 



(T+-T_), T (8) = 2ztt 



cosh7r + T+ - T_T_ 



t 



(3.26) 
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Here the notation was introduced for the Fourier transforms 



^'/A^kf^j. (3,7, 



Using one-loop expression for the coupling constant we get 



t r ^f\ozl/z 2 A 2 \~ 2r±//3 ° 

T ± = -L f d 2 z e"*** , 7 9 / 9 QCD • (3.28) 



r± J ' V logA 2 /A 



2 

QCD 



Let us first consider the properties of the obtained expressions ( |3.26|) and ( |3.27| ) in the leading 
log if: approximation, a s logt/A 2 ~ 1. In this limit we may freeze the argument of the coupling 
constant in ( 3.27|) and perform Fourier transformation to get 



2tt VmV A 2 / r(l - a»r + 



This expression is similar to that (|2.7|) in QED and the only difference is that in QCD the 
s— dependence of the amplitudes is governed by the asymptotic behavior of the eigenvalues 
F±. In the case of abelian gauge group the Wilson lines of figs. 1(a) and 1(b) coincide and the 
cross anomalous dimension in QED is a c-number and not a matrix. As a result, the scattering 
amplitudes have different high energy behavior in QCD and QED. 

Expanding (|3.29|) in powers of a s we find that to the lowest order of PT 



T± = 2a s + 0{a 2 s ) , 

where higher order corrections involve log s and log t terms. After substitution of these values 
into (|3.26| ) we recover the expressions for the scattering amplitude in the Born approximation 

TC°) = + 0(a 2 s ) , = + 0(a 2 ) . 

t mr 

To find the behavior of the scattering amplitude (|3.26|) in the limit of high energies s 3> m 2 , or 
equivalently for large cross angle Q2-5p , we start with the leading logt and logs approximation, 
a s logs/m 2 logt/A 2 ~ 1. We use (3.24) and substitute T+ = N log s/m 2 and T_ = into (|3~29p 
and ( |3.26| ). The result has the standard reggeized form [I8|: 



lffl = 0, TS=T Born .(4) , (3.30) 



m 2 



with the Regge trajectory a{t) = —^Nlog^. We can improve the leading logarithm approxi- 
mation ( |3.30| ) by taking into account nonleading corrections to ( |3.29| ) in the limit a s log s/m 2 ~ 
a s \ogt/\ 2 ~ 1 . Using the large energy behavior of the eigenvalues (3.24) we find the following 
expressions for the amplitudes of singlet and octet exchanges 



T (0) = -T r_(T+ - T_) , T (8) = 2ivrT T+ - — T_ I , (3.31) 




with T = 2^37 and Y± defined in (3.24). The functions T + and T_ have a Regge like behavior 
(|3.29| ) and for T + 3> T_ we have T + C T_. For T_ = the expression (|3.31|) coincides with 
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the leading log result (|3.30| ). We stress that the amplitude of the singlet exchange is 
unvanishing only due to nonleading logarithmic corrections. Moreover, as follows from (|3.31|) , 
the amplitude gets dominant contribution from T_ and not from T + . The same is true for 
the amplitude of the octet exchange in which exponential fall of T + can not compete with 
the factor r_/T + in front of T_. Thus, the high-energy asymptotic behavior of the scattering 
amplitude ( |3.31|) is dominated by the contribution of T_ and as a consequence the amplitude 
of the octet exchange is suppressed by the factor 1/T+ = C(log _1 s/m 2 ) as compared to that 
to the amplitude of the singlet exchange. Comparing this behavior with (|3.30|) we conclude 
that nonleading logarithmic corrections drastically change the leading log asymptotics of the 
scattering amplitude. 

Expression ( |3.29| ) was found under condition that we neglect the r— dependence of the 
coupling constant a s (r) in (|3.27| ). We notice, however, that ( |3.29| ) becomes divergent for the 
special value of the coupling constant a s = 2tt/T±. To understand the origin of this divergence 
we deduce from ( |3.27| ) that it is the distance between quark trajectories, z 2 , which fixes the scale 
of the coupling constant. Then, for the perturbative expansion to be meaninfull this distance 
should be much bigger than the QCD scale 1/Aq CD . On the other hand, to get the scattering 
amplitude in ( |3.27|) we have to integrate in ([3.271) over & U possible z 2 including large distances. 
For the z— integral in ( |3.27| ) to get dominant contribution in PT one has to impose additional 
conditions on T± and g 2 , or equivalently on s and t. Otherwise, the perturbative expression 
(|3.27| ) become divergent at large distances. The analysis of PT applicability conditions for 
( 3.27| ) can be performed similar to the analysis given in [Tj| for the Drell-Yan process. 

The properties of the scattering amplitude are closely related to the particular form of the 
one-loop cross anomalous dimension (3.21). Natural question arises whether the asymptotic 
behavior of the scattering amplitudes will be changed after we take into account higher order 
corrections to the cross anomalous dimension. In the next section we answer this question by 
performing the calculation of the matrix of cross anomalous dimension to two-loop order in the 
Feynman gauge. 



4. Two-loop calculation of the cross anomalous dimen- 
sion 

To find two-loop expression for the cross anomalous dimension we have to calculate the line 
functions W\ and Wi to the second order of PT by taking into account the next terms in the 
expansion of Wilson lines in powers of the coupling constant and effects of gluon self- interaction. 
The resulting Feynman diagrams contributing to W\ are shown in figs. 3-5 and we divide them 
into three different groups 

1. QED-like diagrams of fig. 3 in which gluons don't interact with each other; 

2. Self-energy diagrams of fig. 4 containing self-energy corrections to the gluon propagator; 

3. Three-gluon diagrams of fig. 5 having three-gluon vertex of self-interaction. 

Evaluating the diagrams of fig. 5 containing three gluon vertex, we notice that the Lorentz 
indices of the three gluon vertex T^ up are saturated by the vectors v i or t> 2 depending to which 
line the gluons are attached. As a consequence, the diagrams with all three gluons attached to 
the same line vanish as T ^pV^v^vl = 0. We didn't include these diagrams in fig. 5. 
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There many different Feynman integrals corresponding to the diagrams of figs. 3-5 but, as 
we show in Appendix A, their calculation is reduced to the evaluation of only few "basic" 
diagrams after we take into account the properties of the Feynman integrals. The one-loop 
diagram of fig. 2(1) and two-loop diagrams of fig.3(l)-(5), fig.4(l) and fig.5(l) form the set of 
the basic diagrams. Moreover, among them there are the diagrams which have been already 
calculated in pO . 



4.1. Cusp anomalous dimension 



As was shown in [p0[ , the diagrams of fig.3(3)-(5), fig. 4(1) and fig. 5(1) contribute to the two- 
loop cusp anomalous dimension r cusp (7, g) of Wilson loops. This anomalous dimension appears 
in the RG equation for the Wilson loop if the integration path has a cusp with the angle 7. 

The fact that the same diagrams of fig.3(3)-(5), fig.4(l) and fig. 5(1) contribute to the "cusp" 
and "cross" anomalous dimensions suggests that there might exist a relation between these two 
completely different functions of the angle 7. Indeed, the explicit form of this relation will be 
suggested in sect. 5. Both anomalous dimensions depend on the representation of the SU(N) 
in which the Wilson loop is defined. One might expect that if the relation between them does 
exist it contains r cross and r cusp taken in the same representation. In reality, as we will show in 
sect. 5, the two- loop cross anomalous dimension calculated in the fundamental representation is 
related to the two-loop cusp anomalous dimension defined in the adjoint representation of the 
SU(N). 

The two-loop expression for r cusp (7,g) in the adjoint representation of the SU(N) is given 

i-.v M 



r cU sp(7,s) = -iV( 7 coth 7 - 1) + ( —AM $ cusp ( 7 ) , (4.32) 

where $ C us P (7) = ^qed + $seif + ^3-giuon and different kinds of the diagrams (QED-like of 
fig.3(3)-(5), self-energy of fig. 4(1) and three-gluon of fig.5(l)) contribute to the corresponding 
functions: 

2 P P 
&qed(i) = coth 7 / dtpipyy — ip) cothip — cotli7 / dipip coth ip + 7 coth 7 — ±- 

31 

$seif(7) = ^(7coth 7 -l), (4.33) 
3d 

. . 1 . . , n . P , , Ip COth Ip — 1 7T 2 

$3 g iuon(7) = -q smh(2 7 ) / dip— -2 —2-r ~ 7^j(7coth7 - 1) . 

2 Jo sinh 7 — smh ip 24 

The properties of r cusp (7, g) to higher orders of PT have been studied in [ED]. 



4.2. QED like diagrams 

The same QED-like diagrams of fig. 3 contribute to the line functions W\ and W2 with dif- 
ferent color and combinatorial factors. The explicit expressions for the color factors and the 
combinatorial weights are given in Appendix B. 

To get the contribution of the QED-like diagrams of fig. 3 to the line functions W\ and W2 
we evaluate the corresponding Feynman integrals using the relations found in Appendix A, 
multiply them by color factors and combinatorial weights and sum the resulting expressions. 
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Finally, after lengthy calculation we get the following expression for the contribution of the 
diagrams of fig. 3 to W\ 



QED 



77" 



log 1 



x 2 



N 2 + 1 2 2 Z7T , , ,, . 

-7r coth 7 7coth 7H cotli7 |1) 



S7V 2 



277 

T 



7T 2 , ' ~ ., ,. , ,, / 77 



+ ( -— coth 2 7 + — A^7coth 2 7 - ^-A^ coth 7 J|2) 
47V 8 4 



(4.34) 



We notice that this expression does not contain a single logarithm of \i. Although some in- 
dividual diagrams do contain log(^ 2 /A 2 )— terms they cancel in the total sum of the diagrams. 
In terms of the RG equation ( |2.10| ), this means, that the QED-like diagrams don't contribute 
to the two-loop expressions for the elements Y\\ and Ti2 of the matrix of cross anomalous di- 
mensions. The contribution of the QED-like diagrams to the line function W2 can be found in 
analogous way. The resulting expression has the following structure: 



Aqed log 2 + -Bqed log 



A 



(4.35) 



where in contrast with ll 7 ^ the coefficient -Bqed in front of a single logarithm of \i is different 
from zero. The coefficient Aqed is given by 



A 



QED 



+ 



7T 2 TV 2 + 2 



47V 



+ 



IN 



TV 2 + 1 
87V 2 




-7T 



-27T7 H — 7V 2 7 2 coth 2 7 



N +1 iir + -7V7 ) coth 7 - -Tvl 11) 
47V 8 1 J 4 J 

-7V 2 7 + ^vrj coth 7 + -7V 2 i|2). 



The expression for the coefficient -Bqed is complicated and it contains the integrals Ji, I2 and ^3 
defined in (A.l). We simplify this expression using the identity ( |A.2|) and after some algebra we 
recover the following relation between the coefficient -Bqed and the function $qed contributing 
to the two-loop cusp anomalous dimension (4.33): 



B: 



QED 



-N<t> 



QED 



(7) [|1> - 7V|2>] - ($q ED (7) - ^qed^vt - 7)) [N\l) - |2>] 



(4.36) 



After substitution of the Wilson line (|4.35|) into the RG equation ( |2.10|) we find that the 
coefficient -Bqed contributes to the elements T 2 i and T 2 2 of the matrix of the cross anomalous 
dimension. The identity ([4.36|) implies that the contribution of the QED-like diagrams to two 
different anomalous dimensions, r cross and r cusp , are closely related to each other. 



4.3. Self-energy diagrams 

Calculation of the self-energy diagrams of fig. 4 to the Wilson lines can be easily done using 
the relations found in Appendix A. The total contribution of the diagrams of fig.4 to the line 
function W\ is given by 

WT lf = (^) 2 -coth 7 (^log 2 £ + |log^) (|1> - N\2)) (4.37) 
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and it has the same color structure and the 7— dependence as the one-loop expression ( 3.18Q 
for W\. For the total contribution of the self-energy diagrams to the line function W2 we get 



W, 



self 



5 2 /i 2 31 /i 2 



x 



7V(7Coth7 - 1 - ivrcoth7)|l) - (A^ 2 (7coth7 - 1) - i7rcoth7j|2) . (4.38 



Having the relation (|4.36|) , it is tempting to find the analogous relation between the coefficient 
B Be] f in front of log(/i 2 /A 2 ) in and the function $ se i f defined in (4.33). It can be easily 

checked that the relation does exist and has the following form: 

Sseif = \n<$> BC i f (7) [|1> - N\2)\ - ^($scif(7) - <W^ - 7)) [N\l) - |2)] . (4.39) 



However, in contrast with -Bqed ; the second term in I? se if has additional factor one-half. 



4.4. Three-gluon diagrams 

The diagrams of fig. 5 contribute to the Wilson line W\ with the same color factor defined in 
Appendix B. Then the identity (|A.3|) implies that the total contribution of the diagrams with 
three gluon vertex to the line function W\ vanishes as 

jySgluon ^ F 3gluon + ^Sgluon + ^Sgluon = Q _ 

For the line function W2 the color factors of the diagrams of fig. 5 are different and their total 
contribution is unvanishing 

W 2 3glU ° n = (^) 2 l0g^5 3 gl U on, (4.41) 

where the coefficient B 3g i uon is given by 

£ 3g luo„ = ^3gluon(7) [|1> - AT)2>] - ($ 3g luon(7) ~ $ 3gluo„(^ - 7 )) [N\l) ~ |2>] . (4.42) 

We recognize that this relation between -B3 g i UO n and the contribution of three-gluon diagrams, 
$ 3g i uon , to the cusp anomalous dimension is exactly the same as that for the QED-like diagrams 
in flOSp . 



4.5. Two-loop cross anomalous dimension 



Using the results of calculation of the Feynman diagrams of fig. 3-5 we find the two-loop cor- 
rections to the line functions W± and W% as follows: 



WP = W/ QED + Wf { + H/ 3g 



luon 



1, 2) 



where W® ED was defined in (4.34) and Q43p , Wf { in ( |Q7D and (4.38) and W^ luon in Q440D 
and ([4.41 ). Recall, that the two- loop expressions for the Wilson lines 

(o = l, 2), 



W a = H/ (0) + + H/ (2) 
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should satisfy the RG equation (|2.10|) . Here, is the Born terms defined in ( p,13| ) and 



is the one-loop correction given by ( |3.18|) and (|3.19|) . As a nontrivial test of our calculations we 



find after cumbersome calculation that the two-loop Wilson lines W\ and Wi do obey the RG 
equation ( p. 10 ). Although the coefficients in front of log 2 (/i 2 /A 2 ) in the two-loop expressions 



for the diagrams of fig. 3-5 are complicated, being combined together they give rise to the 
beta-function term in the RG equation ( p.lOj ) for W a . 

The RG equation fl2.10| ) implies that the coefficients in front of log(/x 2 / A 2 ) in the two-loop 

(2) 

expression for the Wilson line Wi contribute to the elements Tn and r\ 2 of the matrix of 

(2) 

cross anomalous dimensions r cross and analogous coefficients in W 2 contribute to the elements 
T 2 i and r 2 2- Notice that Wi and Wi 9luon do not contain log(yU 2 /A 2 ) terms and the elements 
Tn and I\ 2 get contribution only from the self-energy diagrams of fig. 4. Using the expression 
( 4.37Q for Wi we get two-loop corrections to the elements of the matrix of cross anomalous 



dimension 

Ti 2 W n) = - 

36 V 7i 



^nil, g) = -%:(— ) «vrcoth7 



rS 2 2 ) (7^) = §f-V^iVcoth 7 . 



36 V 



7T 



These expressions differ from analogous one-loop expressions (3.21) only by the factor ^-iV||. 

(2} (2) 

Moreover, there is a simple relation between V\i and T^ 2 and the two-loop contribution of the 
self-energy diagrams to the cusp anomalous dimension 

rfiW) = -^(^iv) 2 ($ sclf ( 7 ) - $ sclf (m - 7 )) 



T { ${1,9) = n(^ N ) ( $ sdf(7) - <W^ - 7)) • (4.43) 

(2) (2) 

Calculating two remaining elements, T 21 and T 22 , we find that all three different kinds of the di- 
agrams, W^ ED , W| elf and W2 glwn , contribute to them. Since the corresponding B— coefficients 
are related to the cusp anomalous dimension by the identities ( |4.36|) , ( |4.39|) and ( |4.42| ), the final 



expression for the elements of the matrix of cross anomalous dimension can be represented as 



r 



(2) _J_/\^£ 



21 



rg? - ( 



N V 



~N [$cusp(7) - 2($ cusp (7) - $cusp(w - 7)) + (^solf(7) - ^aclf(«7T - 7))] 
7T / 



*cu S p(7) - ^(*cu S p(7) - $cusp(i7T - 7)) + ^($self(7) - $self(«7T - 7)) 



(4.44) 

where $ C usp = $qed + ^seif + $3giuon wa s defined in (4.33). The terms containing $ se jf appeared 
in the r.h.s. of these relations due to the additional one-half factor in the expression (|4.39| ) 
for -Bscif- The terms containing $ C usp(7) look like two-loop corrections to the cusp anomalous 



dimension Q4.32|) defined in the adjoint representation of the SU(N). 



5. Properties of the cross anomalous dimension 

Having the explicit relations (4.43) and (4.44) between the two- loop corrections to r cross (7,g) 
and r cusp (7,g), we summarize our two- loop calculations by representing the matrix of cross 
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anomalous dimensions in the following form: 

Tcross = -^ r i(7,^) 



1^ = ^(7,2) (5.45) 

r cross = -jj [r C usp(7,^) + Fi(i,g) + r 2 (7,fl')] 

Tcross = r cU sp(7^) + -^J [Fl(j,g) + T 2 (j,g)} , 



where the notation was introduced 



Fi(l,g) 



a s , T 31 f a s „ r 
— N + — — iV 

7T 36 V 7T 



m coth 7 



r 2 (7, a) = -2(r cusp (7, #) - r cusp o - 7, #)) . (5.46) 

We note, that the form of this expression was deduced from the properties of two-loop cor- 
rections (4.43) and (4.44) and it comes as a surprise that the one-loop expression (3.21) for 
r c ross(7,5') also obeys this relation. Interesting property of the obtained expression is that it 
expresses four originally independent elements of the matrix r cross (7, g) in terms of only two 
functions, r cusp and ^(7, g). Written in this form, the relation (5.45) admits natural gener- 
alization to all orders of perturbation theory. Although we do not have a proof that (5.45) is 
valid to all orders of PT, we give below arguments in favour of (5.45). 

Let us examine whether the properties of one-loop r cross found in sect. 2 are still valid after 
we take into account two-loop corrections in (5.45). 

5.1. Imaginary part 

The cusp anomalous dimension T cusp ( / -f,g) takes real positive values for real angles 7 while 
rcusp(^~7 ) g) has a nontrivial imaginary part. It turns out that in the expression for the cross 
anomalous dimension (5.45), 1^1(7, g) and r 2 (7, g) have pure imaginary values and change sign 
under reflection 7 — > m — 7: 

(r a (7,^)) t = -r a (7,#), r a (i7r-7,#) = -T a (j,g), (a = 1,2). 

This property is obvious for the function 1^1(7, g) defined in (5.46). To prove it for the function 
r 2 (7,g), we apply the following property of the cusp anomalous dimension 

2r cusp (7, g) = r cusp (ivr - 7, g) + r cusp (-ivr - 7, g) 
= r cusp (ivr - 7, g) + [r cusp (ivr - 7, g)} ] 

which can be easily verified using (|4.32| ) ,B As a result, the function r 2 (7,g) can be expressed 



as 



r 2 (7, g) = 2i Im r cU s P (wr - 1, g) = r cusp (z7r - 7, g) - r cusp (-27r - 7, g) 

= r cusp (z7r - 7, g) - r cusp (z7r + 7, g) 



1 In fact, the integrals I\, I2 and ^3 defined in (A.l) satisfy analogous relation. 
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where in the last identity we used the fact that r cMsp (7,g) is an even function of the angle 7. 
Let us calculate the determinant of the matrix (5.45): 

N 2 - 1 

detr cross (7,#) = r x (7,(/)(ri(7,</) + r 2 (7,ff)) ■ 

We notice that r cusp does not contribute to this expression due to special form of the matrix 
(5.45). Although the matrix r cross is complex the determinant turns out to be real even function 
of the angle 7, since both Ti and T 2 are imaginary odd functions of 7. 

5.2. Large 7 asymptotics 

Let us consider the asymptotic behavior of the cross anomalous dimension for large angles 7 
between velocities v% and v 2 . We recall that 7 is defined as an angle in Minkowski space- 
time and it may have arbitrary values. The large 7 limit ( |2.5| ) corresponds to the high energy 
scattering of quarks with velocities Vi and v 2 . 

To find the asymptotics of the expression (5.45) for 7 ^> 1 we need to know the analogous 
behavior of the functions r cusp (7, g) and r cusp (i7r — 7, g). To this end we study the most general 
situation and consider r cusp (x, g) with x being the complex quantity with large (negative or 
positive) real part: (Rex) 2 ^ 1- After some algebra we get from the two-loop expression ( |4.32| ) 
for the cusp anomalous dimension that T cnsp (x,g) is a linear function of x in this limit, 

r cusp ( X ,s) = X sign (Re X ) T cusp (g) + 0((Re X )°) , (5.47) 
where the two-loop expression for the coefficient T cusp (g) is given byi 

r -^7" + P) 36 "12 ' 



Moreover, as was shown in [J20| , this asymptotic behavior is valid to all orders of PT. The higher 
order corrections to r cusp (7,g) only modify the coefficient r cusp (g). Using the relations (|5.47| ) 
and (5.46) we get the asymptotics of the function T 2 for large positive angles 7 as 

T 2 ( 1 ,g) = -2tnT cusp (g) + 0( 1 - 1 ) 

and the analogous expression for the function Ti is 

^(7, <?) = mT x {g) + 0{r l ), Ug) = + ^(—nY . 
Thus, the large 7 asymptotic behavior of the matrix of cross anomalous dimensions (5.45) is 

r ii = _j^r ( \ 

cross jy2 ^ ' 

ITT 

rlLs = jjTi(g) (5.48) 

3 %7l 

r cro SS = -7^r cU sp(^) + — (2r cusp (^) - Y^g)) 

r cross = 7r C us P (#) - ^(2r cusp (^) - v 1 {g)) , 



2 We recall, that T cusp is defined in the adjoint representation of the SU (N) group. 
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where we omitted nonleading in 7 terms in real and imaginary parts of the matrix elements. 
Let us calculate the eigenvalues r± of the matrix r cross (7, g). From (5.48) we find their sum 

as 

r+ + r_ = trr cross (7,^) = (7 - -^Jr cross (#) 

and we may expect that both eigenvalues are of order 7. However, calculating the determinant 
of (5.48) we find the product of eigenvalues behaves as 7 , 

2 N 2 - 1 

T+r_ = detr cross (7,#) = tc Ni T l (g)(2T cusp (g) - T^g)) . 

This means that one of the eigenvalues is much larges than the second one: T + = 0(7) and 
T_ = 0( , ~f^ 1 ). The explicit expressions for the eigenvalues are: 



f 2iix 

v 7- ivV 



Vcusp^) (5.49) 



r (i + ^El^ 2 ^ 2 -ir 1 (^)(2r cuap (^)-r 1 (( ? )) 
"'" \i iv 2 7 2 y iv 4 W^) 

where we omitted nonleading in 7 terms in real and imaginary parts of these expressions. 

To get some insight into the general structure of the matrix r cross , let us study the large N 
limit defined as follows: 

a,N = fixed, iV — > 00 . 



We find from (|4.32| ), (5.46) and (5.45) that the anomalous dimensions r cusp (7,g), ^(7, g) and 



1^2(7,(7) survive in this limit. At the same time, all the elements of the matrix r cross vanish 
except of the element 1^2 which becomes equal to the cusp anomalous dimension. As we show 
in Appendix C, this observation is valid to all loop order in PT and the elements of the matrix 
of cross anomalous dimension have the following large N behavior 

r^s = o(n~ 2 ), tZ ss = oiN- 1 ), r 21 oss = o(n- 1 ), r c 2 r 2 oss = r cusp ( 7 , g)+o(N- 2 ) , 

(5.50) 

where r cusp (7,g) is defined in the adjoint representation of the SU(N). We notice that the 
expression (5.45) for the cross anomalous dimension to higher orders of PT is consistent with 
large iV behavior ( |5.50| ). 



Summarizing our two-loop calculations we conclude that the higher order corrections to the 
matrix of anomalous dimensions are organized in such a way, that they preserve the asymptotic 
behavior of the eigenvalues (5.49). This allows us to apply the results of sect. 3. 2 and find the 
higher energy behavior of the scattering amplitudes in QCD. 



6. Conclusions 

In this paper we considered the elastic quark-quark scattering at high energy s and fixed trans- 
ferred momentum t. We have shown that all effects of interaction of incoming quarks with soft 
gluons are factorized into Wilson lines which appear as eikonal phases of the scattering quarks. 
As a consequence, the quark-quark scattering amplitude is represented as an expectation value 
of a Fourier transformed Wilson line evaluated along the integration path which consists of two 
semiclasical quark trajectories separated by the impact parameter in the transverse direction. 
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We evaluated the scattering amplitude using the fact that for zero impact parameter the Wil- 
son line has additional cross singularities. We have shown that the asymptotic behavior of the 
scattering amplitude is governed by the matrix of cross anomalous dimension r cross (7, g) which 
depends on the energy of incoming quarks through the angle 7 between quark velocities. We 
performed two-loop calculation of r cross (7, g) and found remarkable properties of this matrix. 
It turned out that four originally independent elements of r cross (7, g) can be expressed in terms 
of only two functions, the cusp anomalous dimension r cusp (7, g) and 1^1(7, g) which comes from 
self-energy diagrams. Written in this form, the elements of r cross (7, g) admit natural general- 
ization to higher orders of PT. The form of the corresponding expressions was confirmed by 
the large N analysis of the Wilson lines. Higher order corrections are organized in such a way 
that they preserve the asymptotic behavior of the eigenvalues of the matrix of cross anomalous 
dimension. This allowed us to apply the one-loop result for r cross (7,g) to find the high energy 
asymptotic of the quark-quark scattering amplitude. This amplitude can be decomposed into 
singlet and octet invariant amplitudes corresponding to the exchange in the t— channel with 
quantum numbers of vacuum and gluon, respectively. To any order of PT theory the octet 
amplitude gets large double logarithmic corrections while corrections to the singlet amplitude 
are nonleading. However, after resummation to all orders of PT these large corrections are 
summed into the exponentially small Sudakov like exponent with the slope depending on the 
eigenvalues of the matrix of cross anomalous dimensions. As a result, with leading and non- 
leading logarithmic corrections taken into account, the singlet amplitude dominates in the high 
energy asymptotic behavior of the quark-quark scattering amplitude. 
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Appendix A. Basic diagrams 

Let us consider the diagram of fig. 3(1). Using the parameterization of the integration path 
(|3.14j ) we find the following expression for the Feynman integral: 

roo ra± roc pO 

-^1 = iw) ( v i v 2) 2 / doix I d(X2 \ das / c£/3 D(viai — V2^)D{v\a2 — V2OL3) , 

JO JO Jo J-00 

where the D— function was defined in (|3.15|) . It is convenient to introduce the following integrals: 

^1(7) — dip ipcothip 



7 



J 2 ( 7 ) = / dtp tp^-ip) coth ip (A.r 







P , . ip coth ib — 1 sinh 7 
J M = / dip — -2— — — ^-log- 



10 sinh 7 — sinh ip sinh?/' 
Thus defined functions obey the following property 

h{m - 7) + 7/1 (in - 7) = Hi) + (™ ~ 7)^1(7) • ( A - 2 ) 
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Then the final expression for the Feynman integral looks like 



1 fa s 



'47T/X 



2\ 2e 



A 2 

r(i + e)r( e ; 

r(l + 2e) 



r(2e) coth 2 7 



'7(i7r — 7) + 2 r yI\(i'K — 7) — 2(in — 



The cross singularities appear in F\ as a double pole in e. However, the diagram contains a 
divergent subgraph and in order to renormalize F± we have to subtract the cross divergence of 
this subgraph from F±. The corresponding counter term is proportional in the MS— scheme to 
the one-loop expression ( |3.16p 



if E °(7) 



c.-t. 



'47r/i 2 \ £ r(l + e) 



X 2 



As 2 



7(z7r — 7) coth 2 7 . 



To renormalize the contribution of the diagram we add F 1 



QED 



to F^® and subtract poles 



c.-t. 



from the sum in the MS— scheme to find the following renormalized expression 

F QED (7) = (^Wirf* _ 7) lQg 2 ^ + 1 [lh{m _ 7) _ _ 7)/i(7)] log ^!| 

The calculation of the diagram of fig. 3 (2) is analogous to that of F® ED and the final renormalized 
expression for the diagram is 



-,QED, 



a. 



(7)=(^— J coth7|- 7 log — + - [Ji(7) -7] log — 
where the corresponding counter term was taking into account. 



Cusp anomalous dimension 



The diagrams of fig.3(3)-(5), fig.4(l) and fig. 5(1) have been already calculated in [p0| . Using 
the results of the paper p0| we get the following renormalized expressions for the diagram of 
fig-3(3) 



,QED 



f ) h{l) coth 2 7 log^ 



where the integral I2 was defined in (A.l). For the diagram of fig.3(4) we have 

^4 QED = (^) 2 -th 7 I i 7 log 2 ^ + ~ [7 - h (7)] bg ^ } , 

for the diagram of fig. 3 (5) we get 



^QED 



a. 



71 



coth 2 7 



1 2 1 2 A 42 
8 7 bg A* 



^2(7) log V2 



A 2 



for the diagram of fig. 4(1) containing self-energy correction to the gluon propagator 

2/5 ,,2 Q1 ,,2 S 



self 



7T 



48 



, 2 V 

log 



A 2 + 72 



31 , (JL 



log j2 ) 7 coth7. 
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for the diagram of fig. 5(1) 

V , I 

7T 



(7coth7 - 1) + -I 3 (7) sinh(27) 
9o 8 



1 

log^ 



As was shown in [20], all these diagrams contribute to the two-loop cusp anomalous dimension. 



Ward identities 

The Feynman diagrams of figs. 3(7) and 3(4) can be interpreted as propagator and vertex 
corrections to the one-loop diagram of fig. 2(1). As was shown in [20|, there are the Ward 



identities for the Wilson lines, analogous to the Ward identities between Green functions in 
QCD, which relate these two diagrams as follows: 

77QED 77QED 

The same identity is valid for the diagrams of figs. 4 (8) and 4(7): 

ipQED rpQED jpQED 

and for the diagrams of figs. 5(2) and 5(1): 



Ff lf = -^Ff lf (0), 

where the factor one-half takes into account that the self-energy correction is a "square root" of 
the propagator and F* (0) means the value of the Feynman integral F^ elt evaluated at 7 = 0. 

Reflection symmetry: 7 — > m — 7 

The integration paths in figs. 1(a) and 1(b) are defined by two vectors v\ and t> 2 in Minkowski 
space-time. Let us consider the properties of the diagrams under the reflection t> 2 — > — v 2 or 
in terms of the angle between vectors 7 — > m — 7. Under this transformation the integration 
path transforms into itself but the orientation of the t> 2 — line becomes opposite to the original 
one. After inverting of the direction of integration along t> 2 — line as J a dxA(x) = — J b a dxA(x) 
we get formally another diagram multiplied by the factor (— )* with # being the number of 
gluons attached to the t> 2 — line. Thus, taking one of the diagrams of fig. 3, we change the angle 
in the corresponding Feynman integral as 7 — > iir — 7 and multiply the result by appropriate 
(— )* to get the Feynman integral corresponding to another diagram. As a result, we find the 
following useful relations between the Feynman integrals of QED like diagrams 

F QED = F QED {m _ F QED = ^QED^ _ ^ ^QED = ^QED^ _ ^ 

F QED = _ F QED (m _ 7) = F QED (m _ ^ ^QED = ^QED^ _ ^ 

F QED = _ F QED (m _ 7) = _ F QED (m _ ^ ^QED = _ ^ > 

and between self-energy and three gluon diagrams 

F| elf = -Ff lf (i7r - 7), F 2 3gluon = -F? glnon (m - 7) . 

Notice that we don't have the same relations between the color factors of the diagrams, because 
after inverting the direction of the integration along v 2 — line the color indices of the gauge fields 
A(j,(x) = A a (x)t a become anti-path-ordered along the v 2 — line. 
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7^0 limit 

Let us consider the diagram of fig.3(16) in which all gluons are attached to only one of the 
lines. A simple observation is that the Feynman integral of this diagram can be obtained from 
the Feynman integral of the diagram of fig. 3 (3) in the limit V\ — > v 2 or in terms of the angle 
7^0. The same property leads to the following relations between QED like diagrams: 

F QED = F QED (0); F QED = pQED^ ^QED = ^QED^ ^QED = p QEB {Q) 

and between self-energy diagrams 

F self = F sdf( ) ? 

where in the r.h.s. the corresponding Feynman integrals are evaluated at 7 = 0. 



Factorization 

We notice, that in the Feynman diagram of fig.3(20) two gluon propagator are independently 
integrated along the path. As a consequence, the corresponding Feynman integral can be 
factorized into the product of two independent one-loop integrals which have been calculated 
in sect. 2. The same factorization property leads to the following relations: 

pQED _ ^l-loop^2 pQED _ ^pl-loop^2 pQED _ pQED _ pl-loop^l-loop 

pQED _ pQED _ pl-loop pl-loop pQED _ pQED _ pi -loop pi -loop jpQED _ / pl-loop\2 

24 25 2 4 ' 26 27 1 3 28 V 3 / 

Using the expressions for the basic Feynman integrals and applying all these relations we can 
find the Feynman integrals for all the diagrams of figs. 3-5 except of that for the diagram of 
fig.5(3). The calculation of the diagram of fig. 5(3) is based on the following property. It 
turns out that the sum of the Feynman integrals corresponding to the diagrams of figs.5(l)-(3) 
vanishes, that is 

F 3gluon = _ F 3gluon _ ^Sgluon = _ F 3gluon^ + pSgluon^ _ ^ _ ^ ^ 

To show this we notice that three diagrams of fig.5 correspond to all possible positions of two 
gluons on one of the lines provided that the gluons are ordered with respect to each other. 
Moreover, applying the reflection symmetry we may extend the integration over position of 
third gluon from [0, oo) to (— oo, oo) and take one-half of the result. The calculation of the 
Feynman integrals is more simple in the momentum space. Let us denote the momenta of 
gluons as k a , a = 1,2,3 with k\ + hi + k 3 = and k 3 being the momentum of a single gluon 
attached to the v 2 line. Then, taking into account the expression for three-gluon vertex, we get 
the following momentum integral for the sum of the diagrams of fig.5(l)-(3), 

d D h, n .... , (d D k x K 



oc KM - (vi ■ v 2 H) J ■ v 1 )5(h ■ v 2 ) J 



k\ (k. + hnh-v,) 



where two delta functions appear after integration over positions of gluons on v\ and v 2 lines 
from —00 to 00. The momentum integral is proportional to the only vector involved, v f , while 
the prefactor is orthogonal to v±. The product is zero which means that the sum of the Feynman 
integrals corresponding to the diagrams of fig.5(l)-(3) vanishes. 
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Appendix B. Color factors and combinatorial weights 

The relation (|3.20|) implies that the color factors for the diagrams of fig. 3-5 have the following 
general form 

C = dSirfw + CAAi' = Ci|l> + C 2 \2) , (B.l) 

where the coefficients C\ and C2 are different for the same diagram contributing to W\ and 
Wi. Recall, that the relations ( |3.20| ) and ( |B.1| ) are unique properties of the fundamental 
representation of the SU (N) group. 



QED like diagrams 

The general form of the color factors is given by (P-l| ) in the fundamental representation of 
the SU(N). Consider, as an example, calculation of the color factor of the diagram of fig. 3(1). 
This diagram contributes to the both line functions, W± and W2, with the color factors C\ and 
C 2 , respectively, given by 

1 AT 2 — 9 

Gi 2 =(«n^,=-^n)+ i ^i2), 

where we substituted the identity ( |3.20| ) and the states |1) and |2) were defined in ( |3.13| ). 
The calculation of the color factors of the remaining QED like diagrams of fig. 3 is analogous. 
The calculation of the color factors of the QED like diagrams of fig. 3 leads to the following 
expressions 

C\ = C\ = C\ Q = Cl 5 = — 11> - — 12) 

N 2 + 1 1 

C 3 = C n = C 13 = C 20 = 4iy2 |1) - 2jv|2) 

7V 2 — 1 at 2 — 1 

/~t\ _|1\ I _|9\ 

7 — O g — 12 — U 14 — U 22 — C- 23 — C- 24 — C- 25 — |1/ i — \Z) 

N 2 - 1 

°16 ~~ °17 ~~ °18 ~~ u 19 ~~ °26 ~~ °27 ~~ °28 ~~ ^2 I / 

1 /V 2 — 9 

CJ = W|1>, = = — \l) + -^12), 

where the lower index numerates the diagrams of fig. 3 and the upper index is related to the 
Wilson line W\. The color factors of the diagrams contributing to the Wilson line W2 are 

/V 2 — 1 N 2 — '\ 
n 2 — n 2 — n 2 — n 2 — n 2 — n 2 — n 2 — n 2 — zl _io\ 

b 2 — °12 — °13 — °19 — °24 — °25 — °26 — °27 — ^jy I 1 / ^jy 2 l Z / 

TV 2 - 1 

_ _ ^2 _ ^2 _ f~i2 _ ^2 _ ^2 _ ^2 ^-|n\ 

o 3 — o 4 — o 5 — o 6 — o 7 — o 21 — o 22 — o 23 — 4/v 2- ' ' 
C 2 = C 2 = C? 7 = + 12) 
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AN 1 1 4/V 2 ' ' 



1 N 2 4- 1 

^2 _ r 2 _ r 2 _ r 2 ll \ _|_ ___L_ |9\ 

The combinatorial weight of the QED like diagrams are given by 

W 3 = W 5 = W 9 = W n = W 16 = W 18 = 2, W 2Q = W 21 = W 28 = 1 

and the combinatorial weight of the remaining diagrams of fig. 3 is equal to 4. 

Self-energy diagrams 

The color factors of the self-energy diagrams of fig. 4 are defined as follows: 

1 TV N 2 — 1 

c\ = cl = -i|l> + ||2>, c\ = c\ = 

AT 2 — 1 N 1 

C? = W 2 = ^12), C* = Q 2 = f |1> " i|2> 

and the combinatorial weights of the diagrams are 

Wi = W 3 = W4 = 2, U7 2 = 4 . 

Three-gluon diagrams 

The diagrams of fig. 5 contribute to the line function W\ with the color factors 

C\ = C\ = C\ = -2if^(t a t b ) iH t c jfj = — 1|1> + y |2) , 

where / afec are structure constants of the 5*?7 (iV) and the combinatorial weights are 

W\ = w 2 = w 3 = 4. 

For the line function W 2 combinatorial weights are the same but the color factors of the diagrams 
of fig. 5 are different 

N 2 — 1 N 1 

C? = ^yH2>, Cl = -Cl = ^\l)-\\2). 



Appendix C. Large N limit of the cross anomalous di- 
mension 

To find relations ( f5.50| ) we consider the RG equation ( [2.10| ) for the Wilson lines W a and multiply 
the both sides of the equation by 5i>i5j>j. After contacting of the color indices we get from W\ and 
W 2 two Wilson loops which are closed at infinity and have the following large N normalization 

{W 1 t(8 iH 5 rj = N 2 w u {W 2 f(8 iH 5 rj = Nw 2 , (Cl) 
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with w a = O(N ) and a — 1, 2. We notice, that such defined wi is given by the vacuum 
averaged product of two "elementary" Wilson loops calculated along v\ and v 2 lines. The 
simplification of W\ appears in the large N limit after we will take into account the following 
factorization property (vacuum dominance) [[12] 



(O|0iO 2 |O> = (O|0i|O>(O|0 2 |O> x (l + 0{N-' 

valid for two arbitrary gauge invariant operators 0\ and 2 . In our case, 0\ is given by 
Pexp(ig /f^ dsv\A{v ±s) and the expression for 2 is analogous with replacement v± by v 2 . As 
a result, the Wilson loop w\ has a trivial value 

Wl = l + 0(N~ 2 ). (C.2) 

since (0|Ci|0) = (0 1 C?2 10) = 1. Notice, that the first nonleading correction to w\ behaves as 
0{N- 2 ) and not like 0(N^ 1 ). 

After substitution of ( |C.1| ) and flC.2j ) into the RG equation (p . 1 0|) we find the following 
identity in the large N limit 

r^oss + ^l 2 oss w 2 = o(n~ 2 ) . 

The requirement that the both terms in the l.h.s. of this relation have the same large N behavior 
leads together with w 2 = O(N ) to the conditions 

^IL s = o(n-% rZ ss = o(N-i), (c.3) 

which are in accordance with (5.45). Let us consider the projection of the Wilson lines W\ and 
W 2 onto the color structure S^jSj/i 

{Wi)tfS ifj 5 fi = Nw[, {W 2 )tf5 iH 6 fj = N 2 w' 2 , (C.4) 

and repeat the large iV analysis for the Wilson loops w[ and w' 2 - After substitution of ( |U.4j) 
into the RG equation ( |2.10| ) we get 

£S = ^ r cross^i + r cross^2 , (C5) 

where T> denotes the RG operator. The Wilson loop w[ satisfies analogous relation. As in 
the previous case, the Wilson loops w' 2 is factorized in the large iV limit into the product of 
two vacuum averaged Wilson loops evaluated along the paths having a cusp with the same 
angle 7 introduced before. As a result, in the large iV limit w' 2 has cusp singularities which 
are renormalized multiplicatively. This means, that, up to 0(N~ 2 ) corrections, w' 2 should 
obey the homogeneous RG equation with the anomalous dimension ^ L r cusp (7, g) equal to the 
cusp anomalous dimension defined in the fundamental representation of the SU(N) times the 
number of cusps (= 2) with Cf = (N 2 — 1)/ (2N) and Ca = N. Since the Wilson loops w[ and 
w 2 behave as O(N ), this condition together with ( |G.5|) implies that 

rSo, = O(N^), Tl 2 oss = r cusp ( 7 , g) + 0{N~ 2 ) (C.6) 

with r cusp (7,(yf) defined in the adjoint representation of the SU(N). Thus, in the large iV limit 
the matrix r cross (7,g) is given by (|C.3|) and (|C.6|) to all orders of PT. 
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Figures: 




Fig.l: Integration paths (a) and (b) entering into the definition of the Wilson lines W\ and 

W 2 , respectively. 







Fig.2: One-loop diagrams contributing to the line function W\. Solid line represents the 

integration path, dotted lines denote gluons. 
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(1) (2) (3) (4) 



Fig. 4: Two-loop self-energy diagrams. The blob denotes one-loop self-energy corrections to 

the gluon propagator. 
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